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We show how power-law corrections to the Einstein-Hilbert action yield a viable extended theory of 
gravity, passing the Solar-System tests, provided that the power-law exponent n is strictly comprised 
between 2 and 3. We implement this paradigm on a cosmological setting outlining how the main 
phases of the Universe thermal history are properly reproduced. As a result, we find two distinct 
constraints on the characteristic length scale of the model, i.e., a lower bound from the Solar-System 
test and an upper bound found by requiring the existence of the matter-dominated phase of the 
Universe evolution. 

We also show how the extended framework can accommodate the existence of an early de Sitter 
phase. Within the allowed range of characteristic length scales, the relation between the expansion 
rate Hi and the energy scale M of inflation is modified, yielding a value of Hi several orders of mag- 
nitude smaller than the one found in the standard picture, i.e., Hi ~ M^/rUpi. The observational 
implication of this fact is that, quite generally, a tiny value of the tensor-to-scalar ratio r is expected 
in the extended framework, that will go undetected even by future missions focused on cosmic mi- 
crowave background polarization, like CMBPol. The suppression of primordial tensor modes also 
implies that the infiationary scale can be made arbitrarily close to the Planck one without spoiling 
the current limits on r. 

Finally, considering the same modified action, an analysis of the propagation of gravitational 
waves on a Robertson- Walker background is addressed. We find that, in the allowed parameter 
range, the f{R) correction does not significantly affect the standard evolution. 

PACS numbers: 98.80.-k; 98.80.Es;04.50.Kd; 04.25. Nx; 04.30. Tv 
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I. BASIC STATEMENTS 

Several issues in cosmology and quantum field theory seem to suggest that Einstein's General Relativity (GR) 
should be extended in order to cure shortcomings emerging from observations and self-consistent unification theories. 
In early-time cosmology, the presence of the big bang singularity, flatness and horizon problems [l| led to the conclusion 
that Standard Gosmological Model is inadequate to describe the Universe at extreme regimes. On the other hand, 
GR is a classical theory that is not expected to work as a fundamental theory, when a full quantum description of 
spacetime is required. Due to these facts and, first of all, to the lack of a self-consistent Quantum Gravity theory, 
alternative theories of gravity have been pursued in order to find, at least, a consistent semi-classical scheme where 
GR and its successes could be recovered. Extended Theories of Gravity (in particular f{R) gravity) represent a 
fruitful approach and have become a sort of paradigm in the study of gravitational interaction based on corrections 
and extensions of the Einstein scheme Beside fundamental physics motivations, these theories have acquired a 

huge interest in cosmology since they "naturally" exhibit inflationary behaviors able to overcome the shortcomings 
of the Standard Gosmological Model and the associated schemes are, in many cases, able to match the observations 

The appeal of Extended Theories of Gravity also resides in the fact that they could allow to explain the gravitational 
anomalies, usually interpreted in terms of "dark" components of the cosmological fluids, without invoking exotic 
physics, i.e., physics beyond the standard model of particle physics In fact, GR has been successfully directly 
probed only in the weak-field limit (e.g., Solar-System experiments) and also in this case there is room for the 
implementation of alternative theories of gravity which are not at all ruled out, as discussed in many recent studies 
[13.[l3|. In particular, it is possible to show that several f{R) models can satisfy both cosmological and Solar-System 
tests [IJ, ll^ and they can be constrained as the scalar-tensor theories. Moreover, such theories give rise to new 
effects able to explain the so-called Pioneer anomaly (see for example [lB| and references therein) . It is important to 
remark that any modification of the Einstein-Hilbert (EH) Lagrangian is refiected onto a deformed gravitational field 
dynamics at any length scale investigated or observed. Thus, the success of f(R) gravity in the solution of a specific 
problem has to match consistency with observation in different length scales [Ifl]- In this respect, we emphasize that 
the small value of the present curvature of the Universe leads us to believe that, independently of its specific 
functional form, the f{R) term must be regarded as a lower-order expansion in the Ricci scalar. On the other hand, 
it is easily understood that the peculiarities of such an expansion will be extremely sensitive to the morphology of the 
deformed Lagrangian. 

The most immediate generalization of the EH action is to deal with a function of the Ricci scalar R analytical in 
the point i? = 0, so that its Taylor expansion holds 13]. This approach is equivalent to deal with a polynomial form of 
f{R) 0, whose free parameters are available to fit the observed phenomena on different sectors of investigation. 
Despite the appealing profile of such a choice, it is extremely important to observe that it is not the most general 
case, since real (non-integer) exponents of R are, in principle, on the same footing as the simplest case (20l]. 

In this paper, we focus on such an open issue and we will develop a modified theory of the form f{R) — R + qR"^, 
where q and n are two free parameters to be constrained at physical level • In [l^l , it is shown how the details 

of the model lead to deal with rational non-integer numbers for n and to restrict it in the most appropriate interval 
for physical interpretation at low curvature, 2 < n < 3. In fact, by analyzing the weak-field limit of our modified 
Lagrangian and retaining only those non-integer powers in R which represent the dominant effects after the linear 
approximation, the choice 2 < n < 3 allows us to distinguish the dominant non-Einsteinian terms from the non-linear 
ones of GR. The explicit solution of the system derived in [Hj is here re-analyzed confirming the existence of a lower 
bound for the characteristic length scale of the model (defined by means of the parameter q) and of two other typical 
lengths which guarantee that our f{R) model fulfills the Solar-System constraints. 

The cosmological study of the resulting modified Friedmann-Lemaitre- Robertson- Walker (FLRW) dynamics is also 
addressed as the central theme of this paper. The fundamental equation for the scale factor is derived for a perfect 
fluid matter source and the main phases of the Universe evolution, i.e., the radiation- and matter-dominated era, are 
separately discussed. By assuming a generic power-law (in time, t) behavior of the scale factor, we show how the two 
phases are asymptotically properly reproduced. In particular, in the asymptotic limit toward the singularity (taken in 
t = 0), the radiation-dominated scale factor expansion ~ t^^^ is recovered for all values of the FLRW spatial curvature 
and a phase of power-law inflation (^ i"/^) can be generated for vanishing curvature. On the other hand, in the limit 
of large t, the matter-dominated law for the scale factor ~ t^/"^ is reproduced neglecting the spatial curvature. 

The main result of the cosmological analysis is the existence of an upper bound of ^ 70 pc (weakly dependent on 
the value of n for 2 < n < 3) for the characteristic length scale of the model which guarantees the correct matter- 
dominated Universe evolution. Combining this cosmological bound with the Solar-System constraints, we deflne the 
allowed values of the parameter q in our f{R) scheme, thus providing a restricted region of the parameter space in 
which the model can be implemented both in the weak-field limit and in the cosmological framework. 

After the analysis of the power-law Universe evolution, a study of the existence of the standard (exponential) 
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inflationary behavior in our generalized framework is also provided. In particular, we show how an exponential early 
evolution of the Universe is still present in the modified scheme; however, the existence of a minimum length scale 
for the model, as imposed by the Solar-System tests, implies that the expansion rate Hi is much smaller than the 
standard one, i.e.. Hi ~ M^/mpi, by several orders of magnitude. We then discuss the observational implications 
of this modified dynamics, underlying how the production of tensor perturbations is much less efficient than usual. 
This allows to have the infiationary scale M arbitrarily close to the Planck one without spoiling the current limits on 
the tensor-to-scalar ratio r. More generally, we express r as a function of the parameters of the model and of M and 
we find that a tiny value of the tensor-to-scalar ratio is expected in the extended framework. In particular, it is very 
difficult to obtain a value of r > 10~*° unless the exponent n is fine tuned to be very close to 2. Such a small value of 
r will be impossible to detect even for future space-based missions focused on cosmic microwave background (CMB) 
polarization, like CMBPol. 

Finally, we shortly focus on the study of the gravitational wave (GW) propagation on a Robertson- Walker (RW) 
background. The GW equation is derived, using the standard conformal formalism, in the transverse-traceless gauge, 
firstly without specifying the form of f{R)- This way, we obtain a general propagation equation for GWs in the 
case of a flat RW model. We then assume a power-law behavior of the scale factor and our non-analytical form 
of f{R)- We find that the GW propagation is unchanged during the radiation-dominated era, but can be different 
from standard GR, in principle, during the matter-dominated era. However, the existence of a maximum length 
scale implies that the f{R) corrective term can be always neglected during this epoch, and thus does not yield any 
observational consequence. 

The paper is organized as follows. In Section 2, after a brief review of the field equations of f(R) gravity in the metric 
approach, we introduce and discuss our viable f{R) model. In Section 3, we recall and extend the implementation 
of the Solar-System test to our model [13 and derive the lower bound to the length scale associated to the f{R) 
corrections to GR. In Section 4, we carry a study of the resulting modified FLRW dynamics and the inflationary 
paradigm. As expected, this scenario gives us a rather stringent upper constraint to the characteristic length scale of 
the model. In Section 5, we discuss the propagation of cosmological GWs in the presence of f{R) gravity. 

Notations - The signature is set as [ — , -I- ]; Greek indices run form to 3; Latin indices run from 1 to 3; 
(...)' indicates the derivative with respect to R; or (...); denotes the covariant derivative; The dot (...)' is the time 
derivative; (...)_ indicates ordinary differentiation; □ = VV^; We use natural units c = h = 1 and we define x = SttG, 
G being the Newton constant; nipi = 1/VG ~ 10^^ GeV is the Planck mass. 



II. A NON-ANALYTICAL POWER-LAW F(R)-MODEL 

In this paper, we consider the following modified gravitational action in the so-called Jordan frame 

S=^Jd^x^fiR), f{R)=R + qR^, (1) 

where n is a non-integer dimensionless parameter and g < has dimensions [L]^"^^. Such a form of f{R) gives the 
following constraints for n: if i? > 0, all n-values are allowed; if i? < 0, the condition n — £/{2m + 1) must hold 
(where, here and in the following, m and £ denote positive integers). It is straightforward to verify that 5* in Eq.(IT]) is 
non-analytical in i? = for non-integer and rational values of n, i.e., it does not admit Taylor expansion near R — 0. 
Let us now define the characteristic length scale of our model as 

L,{n) ^ |g|i/(2"-2) ^ (2) 

while variations of the total action Stot =5-1- Sm (where Sm denotes the matter term) with respect to the metric 
give, after manipulations and modulo surface terms: 

- i 9^..f - V^V, /' + 5M-n /' - X T^. , (3) 

where T^j^ is the Energy- Momentum Tensor (EMT). 

We can gain further information on the value of n by analyzing the conditions that allow for a consistent weak-field 
stationary limit. Having in mind to investigate the weak field limit of our theory to obtain predictions at Solar-System 
scales, we can decompose the corresponding metric as f/^i, = ry^^ -I- /i^^, where h^i, is a small (for our case, static) 
perturbation of the Minkowskian metric rj^i^. In this limit, the vacuum Einstein equations read 

R,.^- ^V^,^R-miR'''^);^;,^ + nqr]^^nR''-^ =0 , R^3nqnR''-K (5) 

The structure of such field equations leads us to focus our attention on the restricted region of the parameter space 
2 < n < 3. This choice is enforced by the fulfillment of the conditions by which all other terms are negligible with 
respect to the linear and the lowest-order non-Einsteinian ones. 
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III. VIABILITY OF THE MODEL: THE SOLAR-SYSTEM TESTS 



From the analysis of the weak- field limit in the Jordan frame, i.e., Eqs.(l5]), we learn the possibility to find a 
post-Newtonian solution by solving Eqs.dS]) up to the next-to-leading order in h, i.e., np to 0{h"~^), and neglecting 
the 0{h?) contribution only for the cases 2 < n < 3. These considerations motivate the choice we claimed above 
concerning the restriction of the parameter n. The most general spherically-symmetric line element in the weak- field 
limit writes as 



(f -I- <^)dt^ - (1 - *)dr2 - r'^dn^ , 



(6) 



where <i> and 5* are the two generalized gravitational potentials and d^^ is the solid angle element. Within this 
framework, the modified Einstein equations ^ can be integrated [l^ to obtain the following solutions: 



$ = cr + - + $ 
V]/ = * 



n-1 



" I L 



n-1 



A - 

$n 





6nq(3n-4)(n-l) 
(n-2)^ 


2-n 

•} 




6n(3n-4)(n-l)" 


1 

2-n (n-2)^ 




(n-2)2 


6(3n-4)(n-l) 




6n(3n-4)(n-l) 


1 

2-n (n-2) 




[ {n-2Y 


3(3n-4) ' 



(7a) 
(7b) 
(7c) 



where the integration constant 5 has the dimensions [L] and the dimensionless integration constant a can be set equal 
to zero without loss of generality. The integration constant A has the dimensions [L](2»i-2)/(2-n)^ g^^^j (^^^ ^-^^ ^^^.g 
dimensionless, accordingly. Moreover, one can check that and are well-defined only in the case n — {2m + 
while we get A> Q since we assume g < 0. In agreement to the geodesic motion as expanded in the weak field limit, 
the integration constant 5 results equal to 5 ~ , where — 2GM is the Schwarzschild radius of a central object 
of mass M. 

The most suitable arena to evaluate the reliability and the validity range of the weak-field solution ([7]) is the Solar 
System To this end, we can specify Ea. (l7bl) and Eq. ([7cl) for the typical length scales involved in the problem 

and we split 3> and \1/ into two terms, the Newtonian part and a modification, i.e., 

<P = ^N + ^M^ -rfjr + $„(r/L,)2(»-i)/(»-2) , (8a) 
^ = ^N + ^M^ -rg/r + *„(r/L,)2("-i)/("-2) , (8b) 

here, the integration constant S of Eas. (j7b[) - d7cl) is (5 = — Tq = —2GMq {Mq being the Solar mass). While the weak- 
field approximation of the Schwarzschild metric is valid within the range <^ r < oo because it is asymptotically 
flat, the modification terms have the peculiar feature to diverge for r ^ oo. It is therefore necessary to establish a 
validity range, i.e., VMm <C ?" <C VMax, related to n and Lq, where this solution is physically predictive [21']. 

Since we aim to provide a physical picture at least of the planetary region of the Solar System, we are led to require 
that $M and m remain small perturbations with respect to $jv and ^at, so that it is easy to recognize the absence 
of a minimal radius except for the condition r 3> Tq. One can now define a typical distance Lq corresponding to 
the distance where the Newtonian and post- Newtonian terms are of the same order. Thus, for rg ^ r <C Lq, the 
system obeys thus Newtonian physics and experiences the post-Newtonian term as a correction. Another maximum 
distance Lq* can be futhermore considered, following from the request that the weak-field expansion © should hold, 
regardless of the ratios '^m/'^n and '^m/'^n- Lq and Lq* result to be defined as functions of n and Lg and they 
write as [Tl[2^ 

n— 2 2n — 2 n — 2 

L*Q ^ |rg/$„|3;^L,3"-4 , L*Q* ^ Lg/\<i>r.\2^ , (9) 

and, for the validity of our scheme, Lq Tq must hold, i.e., Lq >• rQ|$„|("~^^/(^"^^^ 

Neglecting the lower-order effects concerning the eccentricity of the planetary orbit, we can deal with the simple 
model of a planet moving on circular orbit around the Sun with an orbital period T given by T = 27r(r/a)^/^, 
where a = {d^/dr)/2 is the centripetal acceleration. We can now compare the correction to the Keplerian period 
Tk = 2'kt'^I'^{GMq)~^I'^ , with the experimental data of the period Te^^p and its uncertainty bT^xp. We then impose that 
the correction arising from the $m term in Ea. (l7bl) is smaller than the experimental uncertainty, i.e., 6Tf,xp/Tf,^p ^ 
\Tk — T\/Tk. Specifying now our analysis for the case of the Earth [8], i.e., T^xp — 365.2563 days and ST^xp — 
5.0 • 10"^" days, we can get a lower bound Lq > Lq^™ for the characteristic length scale j22], as function of n. For 
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our model, from Eas. (j7bp . it writes 



- Mi: 



in) 



1.3689 X 10' 



|$„| 71-1 



n-2 



3n-4, 



(10) 



pc is the mean orbital distance of Earth from the Sun. The minimum length L^^™ mcreases 

2 it tends to 

-^Mtn _ io-6(y^ _ 2). For a typical value n ~ 2.66, one gets - 4 x IQ-^ pc. We 



where rp ~ 4.8482 x 10 

monotonically for 2 < 7i < 3; it tends to a finite value L*^™ ~ 0.05 pc when n ^ 3, while for n 
with the asymptotic behaviour L*^™ « 10~^(n — 2). For a typical value n ~ 2.66, one gets L^*" ^ 4 x lO"'^ 
remark that L^™, by virtue of Eg. dTbl) . is defined only for n — {2m + l)/£. 

Our analysis clarifies how the predictions of the corresponding equations for the weak-field limit appear viable in 
view of the constraints arising from the Solar-System physics. Indeed, the lower bound for Lg does not represent a 
serious shortcoming of the model, as we are going to discuss in Section UV CI 



IV. COSMOLOGICAL IMPLEMENTATION OF THE MODEL 



In order to study how our f(R) model affects the cosmological evolution, we start from the modified gravitational 
action ([T]) and we assume the standard RW line element in the synchronous reference system, i.e., 

^ds^ = - a{tf [ drV(l - Xr2) + r^df^^ ] , (11) 

where a(t) is the scale factor and K the spatial curvature constant. With this form of the metric, the 00-component 
of Eq.Q results, using the Bianchi identity, to be the only independent one and writes as 

/'i?oo + i / - 3(d/a) /" i? = -X Too . (12) 

We assume as matter source a perfect-fluid EMT, i.e., T^i, ~ {p + p) u^u,^+pg^i,, in a comoving reference system (thus 
Too — p)i where p is the thermostatic pressure, p the energy density and denotes the 4- velocity. Assuming the 
equation of state (EoS) p = wp, the p = component of the conservation law T^. = gives the following expression 

for the energy density: p = Po [a/ao]~^(^+™^. 

Using now f = R + qR" with q < 0, we are able to explicitly write Eg. (fT2|) : 



2x0^-3"' -I- G^-nqa^^^^'di-K -d^ -ady-^ + 

+ [ - 6X - 6 d^ + 6" g a^fi-") {-K - d^ - a d)"] + (13) 
+ 6'\n- l)nqd a^^^-") (-K - d^ - a d)""^ [ - 2 d^ - 2X d + a d d a] = , 

where x = XPa'^a^^^^ ■ Let us now assume a power-law a = [i/^o]^ for the scale factor and, for the sake of simplicity, 
we set d = a^t^^^ (clearly, [d] = [L^"^]). Here and in the following, we use the subscript (...)o to denote quantities 
measured today. In this case, Eg. dTSt can be recast in the form 



-6d2i^i2._6dVi4"-2^gd4f4- (c'it-2_6d-2xi-2-)"+ 2xdi-3-t-(i-3-) = 

where Ci = 6x(l - 2a;), C2 = [2;(2?i - 1) - 1], C3 = xd^ix + 2n - 3)(2x - 1), and 
Ci =xd^{2x- 1). 



A. Radiation-dominated Universe 



Here, we assume the radiation-dominated Universe EoS w — 1/3 {p ^ a^^). In the following, we will separately 
discuss the three distinct regimes for x < 1, x > 1 and a; = 1 in the asymptotic limit t ~> 0,. 

In the case x < 1, all terms containing explicitly the curvature K of Eq.(fT4|) results to be negligible for t — >■ and 
asymptotic solutions are allowed if and only if x ^ n/2 which, in the case 2 < n < 3 we are considering, is always 
satisfied. The leading-order term of Eq. (|T4| writes as 



qd^C^ [1 - {C3/Cl)nx^d^] t^^-^n ^ , 



(15) 
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and x = 1/2 and x ~ [2 + 3n — 2n^ ± {4 + 8n + — 12n'^ + 4n^)^/^]/2n are the solutions. Such second expression 
results to be negative or imaginary for 2 < n < 3 and must be excluded. Thus, the only solution for a; < 1, in the 
asymptotic limit for t 0, is the well-known radiation dominated behavior a ^ t^/^. In the other two cases, i.e., 
for a; 5^ 1, it is easy to recognize that no asymptotic solutions are allowed. Therefore, the approach to the initial 
singularity is not characterized by power-law inflation behavior when spatial curvature is non- vanishing. 

Let us now assume a vanishing spatial curvature in Eg. ([T4| . In can be show how, for K = 0, the radiation- 
dominated solution with w = 1/3 and x = 1/2 is an exact solution (non-asymptotic and allowed for all n- values) with 



Po = 3/{4:xtg), thus matching the standard FLRW case. In the case x > 1, the leading-order terms of Eg. p^ read, 



for t and K = 0, 



qa^C^ [ 1 - {C3/Ct)nx'a' ] t*""^" + 2x = . (16) 

Three distinct regimes have to be now separately discussed. For x > n/2, the leading order of the equation above 
does not admit solutions since it writes simply 2% = and, for x < n/2, the solutions of Ea. (|16|) are those obtained 
in the case for a; < 1. Instead, for x = n/2, and defining Ho — {n/2)/to, one gets 



p, = , p„{n) - ^ (1 - n)("-i)n"(n/2)2-2" [4n -|- (6 - 5n)n' - 4] , (17) 



where we have introduced the dimensionless parameter = Hg "~^ q. We remark that the constraint n = (2™ -|- 
l)/(2^-|-l) (which is in agreement with respect to the one obtained from Solar-System test) must hold in order to have 
Po > since we have assumed g < and therefore go < 0. The function po results to increase as n goes from 2 to 3 and, 
in particular, one can get 216 < Po < 21 024. Finally, for a: = 1 and K = 0, Eq.l^^ reads [1 - n{2n - 2)] t^-'^"' = 0, 
giving n = [1 ± VS\/2. As the previous case, the regime a; = 1 does not admit solutions in the region 2 < n < 3. 



B. Matter-dominated Universe 



Let us now study the matter-dominated Universe EoS w — {p ^ a~^). As previously done, we analyze the three 
distinct regimes for x < 1, x > 1 and x = 1, and, in the limit for t — >■ oo, it is easy to recognize that there are no 
power-law solutions in all these cases for K ^ 0. Setting K = 0, the x ^ 1 regimes do not provide any power-law 
form for cosmological dynamics either. On the other hand, for x < 1 and assuming zero spatial curvature in Eg. (1141) . 
we get the following equation: 

[-6x^a^] + 2xd = a^qC^[-l + Cza^nx' jCW i*"-^" _ (18) 

Since 4x — 2 > 4x — 2n, the term on the right hand side can be neglected in the limit of large t and the equation above 
admits three distinct situations: x < 2/3, x > 2/3 and x = 2/3. Both cases with x 7^ 2/3 do not admit solution. The 
case X = 2/3 admits instead an asymptotic solution for t — 00. In fact, Eq. p8)) reduces to 8a^ = 6% and the FLWR 
matter-dominated power-law solution a — at^l"^ is reached setting po = 4/(3xio). In conclusion, we can infer that, for 
f{R) = R + qR^, the standard matter-dominated FLRW behavior of the scale factor a ^ t^l^ is the only asymptotic 
(as t — >■ c») power-law solution. 

As shown above, the matter dominated solution a ^ t^l^ is obtained ioY K — ^ and asimptotically for large t. In 
order to neglect all the ii'-terms in our /(i?) model, we start directly from the expression of the Ricci scalar [23l |. 
Using a power- law scale factor, we get the t-range (if x 7^ 1/2 and x < 1) 

t « I [x(2x - l)]/[X/a2] I . (19) 

For the matter-dominated era and using standard cosmological parameters 17], one can get the upper limit 
i4r/a2< 0.006 (i7(,)^/^. Thus, setting x = 2/3, we get the bound t < 235/i/o, independently of the form of 

At the same time, if we set x = 2/3, the asymptotic solution po — 4/(3x^o) reached neglecting the right hand 
side (<C 1) of Eq. p^ . i.e, if t is constrained by the following lower limit (wc remind that g^, = q) 

i » Kn,qo)lHo , ^i{n,qo) = \ g^ [ - (4/3)" + 2(2"+i)3-" n(2n - 7/3)] |'/'^""'^ . (20) 

Let us now recall that the matter-dominated era began, assuming Ho^ — 4.3 x lO^^s, at tEq — 5.1 x lO^^/iJo- In this 
sense, we can safely assume p(n,qo) ^ 5.1 x 10~*, which implies an upper limit for |go|, i.e., \qo\ ^ |go|*^°^, where 

|g„|A^"^(n)= [5.1 X 10-*]^^^'"^ I - (4/3)"-f 2(2"+i)3-"n(2n-7/3)P^ . (21) 

It is easy to check that the function |go|*^°^(n) is decreasing as n goes from 2 to 3, in particular, one gets: 
10"^^<|go|*'^'"'< 10"^^ 
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C. Characteristic length scales of the model 



As already discussed in Section 2, the parameter q has dimension [L]^"^^. We have therefore defined a characteristic 
length scale of the model as Lq(n) = jgl^/^^""^). Assuming f{R) corrections to be smaller than the experimental 
uncertainty of the orbital period of the Earth around the Sun, the lower bound PU)) for Lq{n) was found. In order to 
identify the allowed scales for our model and in view of the upper constraint on the parameter = iJ^ q derived 
in the cosmological framework, we can now define the upper limit for Lq{n) as 

if (") = [k„|^"'^1^/('"-'Vir„ , (22) 

which, considering Eq. (PT|) . yields to the constraints 65.59 pc < L^^'^'-^ < 78.37 pc, in the case 2 < n < 3. Assuming 
~ 4.2 X 10^ pc, the two bounds for the characteristic length scales here discussed, i.e, Eq.([TU]) and Eg. ([2^ . are 
plotted in Figure [T] (A). At the same time two other typical lengths have been outlined in Eq.(|9l) for the Solar System. 




Figure 1. Panel A: L^^*" of Ea. (fTO|) and L,^"^ of Ea. (l22l) . The gray zone represents the allowed characteristic-length scales 
of the model. We stress that L^q" is defined only if n = (2m + V)/l, as represented by the dotted line. Panel B: Lq and L*q 
of Eq.([9]). The gray zone represents here the allowed distances for the model. 



Lq represents the minimum distance to have post-Newtonian and Newtonian terms of the same order, while Lq* was 
defined according to the request that the weak-field expansion holds. Setting now n = 23/9 ~ 2.55, one can show 
from Eg. p^ and Eg. ip^ that the allowed scales are 0.0013 pc<Lg< 71.72 pc. In this range, Lq and Lq can be 
plotted as in Figure [1] (B). 

Summarizing, our analysis states a precise range of validity for the power-law /(i?) model we are considering. 
Indeed, for a generic value of n {i.e., not close to 2 or 3) the fundamental length of the model is constrained to 
range from the super Solar-System scale up to a sub-galactic one. Therefore, in agreement with Ea. dSal) . we have to 
search significant modification for the Newton law in gravitational system lying in this interval of length scales, like 
for instance, stellar clusters (see also [23|). 



D. The inflationary paradigm 

After discussing the power-law evolution of the Universe proper of the radiation- and matter-dominated eras, we 
now analyze the inflationary behavior characterizing the very early dynamics. In this analysis, we assume that the 
inflationary expansion is driven by some component {e.g., a scalar fleld) with negative pressure p = —p, as in standard 
inflationary models. In other words, we are interested in studying how the modification of the GR action ([T]) affects the 
de Sitter dynamics (and, indeed, if such a phase can exist in the extended framework). For an interesting approach, 
where the infiationary scenario is instead directly embedded within the modified gravity scheme, see [25l. |26I|. In this 
respect, we hypothesize an exponential behavior for the scale factor of the Universe where 
Hi > and a — aie~^'*^' (here, ti denotes the end of inflation). In the following, we concentrate the attention on 
the solution for vanishing spatial curvature K = Q and, in this case, Eg. (jl3|) rewrites as 

e^"''[q{-l2)'' Hf\l - n/2) - 6Hj] + 2x(a e^^*)!-^- = q . (23) 

Assuming now w = —1, i.e., p — pi ~ const., the equation above reduces to 

[(-1)"12" q{l - n/2)]Hf'' - 6H] + 2xpi = . (24) 
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During inflation, pi can be assumed of the order of the Grand Unification energy-scale, z.e, pi ~ (lO""^^ GeV)^. The 
equation above yields an algebraic relation between Hj (the Hubble constant during the de Sitter phase) and pi , exactly 
as in the standard case; in fact, the latter is recovered for q — and the well-known relation Hi = Hj = y/xPiJ^ 
holds. 

In order to integrate Eq. (p4)) . we focus on a particular value of the power-law f{R) exponent, e.g., n — 29/13 ^ 2.23. 
Using Eq.dini), we obtain L^'" ~ 1.5 x 10"^ pc and, having in mind that Lg = |g|i/(2n-2) .^^g |g| > |^|M»n ^ 

1.4 X 10~^^pc^^/^'^. Let us now fix the parameter q to a reasonable value like q* ^ — 10'^|g|^™ (we remark that 
such assumption will be physically motivated in the next Section). In this case, the solution of Eg. ([Ml) is Hj ~ 
1.8 X 10^^ pc~^ ~ 1.2 X 10~^° GeV. It is worth noting that, for typical values of q like the one used above, the equation 
above admits solutions only in the case n = {2£ + l)/(2m -I- 1), consistently with the previous analyses. 

This study shows how an exponential early expansion of the Universe can still be associated to a constant vacuum 
energy, even for the modified dynamics. However, we see that the rate of expansion is significantly lower than the 
standard one by more than 20 orders of magnitude. In fact, the standard value is Hj = y/xPi/^ — 3.7 x lO^^pc"^ ~ 
2.4 X lO^'^ GeV. Although our estimation relies on the assumption that pj is of the order of the Grand Unification 
energy-scale, nevertheless the values of Hj remains many order of magnitude below the standard value even when 
changing pj, as long as this assumes physically motivated values, as shown in Figure [21 Assuming that the minimum 
number of e-folds required to solve the standard model paradoxes is the usual one, this implies that the duration of 
the de Sitter phase should be much longer (time- wise). Despite this difference, we have shown that it is still possible 
to arrange the cosmological parameters in order to have a satisfactory inflationary scenario. 




(B) 



~ 2x10"" 
I L, ~ 6x10"" pc 



JQ-20 JQ-15 JQ-IO JQU JQ5 JQlU JQ 



10 ,nl5 



p]" IGeV] 



p]" IGeV] 



Figure 2. Panel A; Numerical solution of the modified Friedmann equation (|24|l setting n = 29/13 and Lq = 2.5 x 10 ^ pc, 
corresponding to g = — 10''|g|^^'" (solid line). We also show the standard solution Hi = \fxinj?> (dashed line) and the solution 
obtained assuming that the /(-R) correction term is the leading one (dotted line), for a discussion see Eq. (|27I) . Panel B: The 
same as the panel A, but for (n — 2) ~ 2 x 10~® and ~ 6 x lO^^^pc (again corresponding to g = — 10^|g|*^™). 



Let us now analyze in more detail the implications of the modified Friedmann equation (|24p . First of all, we assume 
that the first term in this equation dominates over the second one, i.e.. 



Hi > (n- 2) ^(^L-i 



(25) 



neglecting a numerical factor of order unity. We will check a posteriori the parameter range for which such approxi- 
mation is reasonable. With a little algebra, Ea. (lM[) can then be recast as 



(-12)"-^ \q\ (n-2)i7f -XP//3. 



(26) 



The positivitcness of the right-hand side requires that (—1)" ^ = 1 and then that n is of the form {2£+ l)/(2m-|- 1), 
as anticipated above. In this case, the relation ([26[) can be explictly inverted to give: 



Hr 



12"|g|(n-2) 



12" (n 



(27) 



Denoting with rriq = the energy scale associated to the parameter q and with AI = p'J'^ the energy scale of 
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inflation, and recalling that x = ^Tr/mpj we can write (neglecting a numerical factor of order unity): 



Hr 



This should be compared with the standard value H j 
velocity thus roughly amounts to 



n-l 



(28) 



\fxpiji ^ M^/mpi. The suppression factor for the expansion 



Hi 



(n-2)- 



AP J 



For consistency, we should check that the solution (|27|) satisfies the condition ([25]) . This is the case if 

pj = > 



1 



327r(n - 2) 



L i 



(29) 



(30) 



We are going to argue that this is practically always the case, due to the existence of the lower bound (|T0|) . In fact, 
since Lq > L^™(n), one can argue that, if the above condition is satisfied for Lq — Lq^"^{n), then it is automatically 
satisfied for all the values of Lq allowed by the Solar-Systems tests. We have said that this is practically always 
the case, because the right-hand side of Eq. (j30p blows up when n — >■ 2, thus in this limit it is increasingly difficult 
to satisfy the condition on the energy density expressed by this equation. This is related to the fact that the f{R) 
correction to the Friedmann equation (1241) during inflation vanishes for n = 2. However, we will show in a moment 
that, for reasonable values of p/, this is not a problem unless the value of (n — 2) is unnaturally tuned to be very small. 
When n is not too close to 2 (for definiteness, let us take {n~ 2)> 10^^) the right hand side of Ea.(pO|. calculated 
at Lq = L*^™(n), lies roughly between (IQ-'^GeV)^, for n = 3, and - (lO^^GeV)^, for n = 2.01. The energy scale 
of inflation is poorly constrained by observations; however, inflation certainly has to involve some physics beyond the 
standard model, and we can safely take M > 1 TeV. By doing this, we see that the condition (PH)) is always satisfied 
as long as n is not very close to 2. More quantitatively, we find that (n — 2) > 2 x 10^"^ is enough to ensure that 
M > 10 MeV. Although [n — 2) could be very small in principle, nevertheless having (n — 2) ~ 10^"^ or smaller seems 
like an undesirable fine tuning of the model. It can be easily understood that the amount of fine tuning required in 
order to spoil the condition (p0| gets larger when considering higher values of the infiation energy scale. The above 
reasoning is confirmed by the inspection of Figure |31 where we show the allowed values of M, the energy scale of 
inflation, as a function of n. Finally, we also note that the condition (|30l) implies, through Eq. ([29| . that in this regime 
Hi /Hi < 1. 




Figure 3. Panel A: In grey: region in the M — n plane where the /(i?)-term in the modified Friedmann equation dominates 
during inflation, irregardless of the value of Lq (as long as Solar-System bounds are satisfled). We remark that this term 
could dominate also in the region below the curve, depending on Lq. Panel B: The same as panel A, but this time in the 
M — log(n — 2) plane, in order to highlight the behaviour for n — t- 2. 



We conclude this section by discussing an important observational implication of the smallness of Hi during inflation 
with respect to its GR value Hi . It is well known that tensor perturbations hij to the background metric are produced 
during infiation. These perturbations physically correspond to GW fiuctuations and give origin to a cosmological 
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background that could be, in principle, detected through measurements of the CMB temperature and polarization 
spectra. The amplitude of the tensor perturbations is expressed through their power spectrum Py(fc), quantifying the 
variance of fluctuations with wave number k. In standard GR, this is found to be 



Prik) = 4x ( ^ 



(31) 



k=aH 



where the k dependence comes from the fact that the right-hand side is evaluated at the time the mode of interest 
leaves the horizon (fc = aH). Since H is nearly constant during inflation, the value of H at horizon crossing is, in a 
first approximation, the same for all fc-modes; this implies that the spectrum should be nearly scale- invariant. In the 
following we will neglect the scale dependence of the spectrum and always take H = Hj. It has been shown that, 
in the framework of modified theories of gravity, the amplitude of the spectrum is modified as follows: 



Pt = 



Ax 



irm 



Hi 

27r 



(32) 



The f'{R) term is given by: 



f'iR) = 1 



nqR""-^ = 1 



(33) 



where we have used the fact that during the de Sitter phase R — 12Hf, and we have substituted expression ^7} 
for Hj. It is straightforward to verify that the condition pop ensures that the second term in the equation above is 
always much larger than unity; thus we have 



np) 



327r 



> 1 



(34) 



This implies that the amplitude of the primordial spectrum of tensor modes is strongly suppressed, first because of 
the large factor f'{R) appearing in Eq. ([32l) . secondly because the value of Hi is itself suppressed with respect to its 
standard value, as discussed above. Substituting the expressions ([27]) and ([34|) for /' and Hi into Eq. ([32l) yields: 



12n7r2 L2 




(n ~ 2)m: 



9 ? 



327rM4 



(35) 



It is interesting to note that, contrarily to the usual, the amplitude of tensor modes is smaller for large values of M , 
and that in the limit n — > 2, it does not depend on M but only on m^, i.e., PT\n-^2 — (l/67r)(m2/mpj). It is useful to 
write the tensor spectrum in terms of the prediction of standard GR, given by Eq. (j31[) evaluted for Hi = Hi, times 
a factor /3, that will turn out to be very small. By doing this, we obtain 



Pt = 



(n-2) 



327rAf4 



Pi 



PPi 



(36) 



where /3 = (3{n, ruq, M). The amplitude of the spectrum is clearly maximized in the case rriq — 1/L^ (n), i.e., when 
the Solar-System bounds are saturated. As discussed above, this maximum allowed mass is a decreasing function 
of n, formally diverging for n 2, but still cannot be made very large unless n is fine-tuned to be very close to 2. 
For n — 2.01, we get the tiny value L^^™ = (5 x lO"^^ GeV)""'^, and it is easy to check that, in order for the factor 
mpiniq/M^ , that roughly controls the magnitude of /3, to be of order unity, one would need M to lie at or below the 
MeV scale, that is completely unrealistic. This is made clearer in Fig. |4l where we show curves of constant (3 in the 
{M, n) and (Af, log(n — 2)) planes. We can conclude that reasonable values of n and M yield Pt ^ Pt- 



The amplitude of tensor modes is measured by the ratio r of tensor perturbations to scalar (density) perturbations, 
evaluated at the pivot wavenumber k^ = 0.002 Mpc~^, i.e., r = {Pt / Ps)\k=ka- The primordial amplitude of scalar 
perturbations is well measured by CMB experiments to be Ps|fe=fco — 2.4 x 10~^; using this, Eg. ips]) can be readily 
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2.0 2.2 2.4 2.6 2.8 3.0 IQ-" lO"' 10"^ 10"' 10"- 10"' 10" 

n n—2 



Figure 4. Panel A: Curves of constant /3 in the (M, n) plane. The value of /? decreases going from left to right and from 
bottom to top. From left to right, the plotted curves correspond to /3 = {10"^°, 10"^°, . . . , 10"^^°}. In aU cases ruq is assumed 
to saturate the Solar-System bounds, in order to give the largest possible value for /3. Panel B: The same as panel A, but in 
the (M, log(n - 2)) plane. From left to right, the plotted curves correspond to /3 = {10"^", 10^^°, . . . , 10"^°°}. 

turned into an expression for r. The tensor-to-scalar ratio is constrained by the recent WMAP data to be r < 0.36 
at 95% confidence level (c.L); inclusion of other datasets lowers the bound down to r < 0.20 at 95% c.l. [H, [2^. The 
present constraints are essentially driven by the temperature and temperature-polarization cross-correlation data [30| . 
The incoming data from the Planck satellite will likely improve the 95% c.l. bound to r < 0.04 [3l|. Tighter constraints 
could be obtained through better measurements (meaning either a detection or a strong upper bound) of the so-called 
B-modes of the CMB polarization: in principle, a future space-based mission focused on CMB polarization, like 
the CMBPol mission currently under study, could detect tensor modes at 99% c.l. for r>0.01 or provide an upper 
limit r< 0.002 at the same c.l. 32j. However, expression (I55t yields a value of r ~ 10~^° for n — 2.01 (practically 
independent from M), and still lower values for larger n (see Figure [5]). This implies that the extended theory of 
gravity considered here predicts a cosmological background of gravitational waves that will be beyond experimental 
reach even in the far future, at least for the very large wavelengths probed by CMB experiments. This is made even 
more clear by the plots in Figure [SI where it is evident that "large" {i.e., potentially detectable) values of r are only 
possible in a tiny region of parameter space corresponding to unnaturally small values of n — 2. More precisely, since 
in the limit of small n — 2 the asymptotic behaviour of r is r ~ 10~^'^(n — 2)~^ (independently of M, and optimistically 
assuming that Lq — L*^*"), it is straightforward to see that, in order to have a tensor-to-scalar ratio at a level of lO"'^ 
or larger, n — 2 has to be of the order of lO"**" or smaller. 



The smallness of the tensor-to-scalar ratio can be considered as a well-definite, falsifiable prediction of the model. 
The immediate implication of the discussion above is that a detection of a non-zero value of the tensor-to-scalar ratio 
in near future would rule out the model considered here, unless one is willing to admit that n is fine-tuned to be 
incredibly close to 2, but yet, not exactly equal to 2, as this value is not admitted by the theory. 

V. PROPAGATION OF GRAVITATIONAL WAVES ON A RW BACKGROUND 

The design and the construction of a number of sensitive detectors for GWs is underway today. There are some laser 
interferometers like the VIRGO detector, built in Cascina (near Pisa, Italy) by a joint Italian- French collaboration, the 
GEO 600 detector built in Hanover (Germany), by a joint Anglo-German collaboration, the two LIGO detectors built 
in the United States one in Hanford (Washington) and the other in Livingston (Louisiana) by a joint Caltech-MIT 
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Figure 5. Tensor-to-scalar ratio r as a function of the inflationary energy scale M, for different values of the index n. Solid line 
n = 2.01; dashed line n = 2.5; dotted line n — 3. In all cases is assumed to saturate the Solar-System bounds, in order to 
give the largest possible value for r. 
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Figure 6. Panel A: Curves of constant r in the (M, n) plane. The value of r decreases going from left to right and from 
bottom to top. From left to right, the plotted curves correspond to r = {10"*°, 10"*^ . . . , 10"^^^}. In all cases iriq is assumed 
to saturate the Solar-System bounds, in order to give the largest possible value for r. Panel B: The same as panel A, but in 
the (M, log(n - 2)) plane. From left to right, the plotted curves correspond to r = {10^™, 10"^^, . . . , 10"^"°}. 



collaboration, and the TAMA 300 detector, in Tokyo (Japan). Many bar detectors are currently in operation too, 
and several interferometers and bars are in a phase of planning and proposal stages (for the current status of GW 
experiments see [HI). The results of these detectors will have a fundamental impact on astrophysics and gravitation 
physics. There will be a huge amount of experimental data to be analyzed, and theorists would face new aspects 
of physics from such a data stream. Furthermore, GW detectors will be of fundamental importance to probe the 
GR or every alternative theory of gravitation [s^ . A possible target of these experiments is the so-called stochastic 
background of GWs. The production of the primordial part of this stochastic background (relic GWs) is well known 
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in literature starting from [35|, in which, using the so-called adiabatically-amplified zero-point fluctuations process, 
it has been shown in two different ways how the inflationary scenario for the early universe can, in principle, provide 
the signature for the spectrum of relic GWs. 

The GW equations in the transverse-traceless gauge write D/i^ = and they are derived from the Einstein field 
equations deduced from the EH Lagrangian C — R. Clearly, matter perturbations do not appear in the GW equations 
since scalar and vector perturbations do not couple with tensor perturbations in the Einstein equations j2j. For a 
general derivation of the propagation equation of GWs in the extended framework, we start considering Eq. ([1} without 
specifying, in the first instance, the form of f{R). For the sake of simplicity, we discard matter contributions and a 
conformal analysis will help to this goal. In fact, assuming the conformal transformation 

gM-=e2^.9M- with e^^ = riR), (37) 
where (j) is the conformal scalar field, we obtain the conformally equivalent HE action 

S = ^Jd^x^[R + C{<P,cP,^)\ . (38) 

Here C (</>, 0;^) is the conformal scalar field contribution derived from 

R^.. - + 2((/.;^0;, - 5m.</);5</)^* - - ;5) , R^ c'^^ (R - 6n0 - 6</);5</)'*) . (39) 

Starting from the action (j38|) and deriving the Einstein-like conformal equations, expressed in the conformal metric 
g^ivi the GW equations writes as 

□ /i^ = . (40) 

Since no scalar perturbation couples to the tensor part of GWs, one can easily recognize that the relation h\ = 
9^^ = g^^ e^'^Sgu = hi holds, which means that hi is a conformal invariant. As a consequence, the plane-wave 
amplitude 

hi = hit) 4 e*-' , (41) 
where el is the polarization tensor, are the same in both metrics. In any case, the d'Alembert operator transforms as 

□ = (□ + 20^^Va) , (42) 

and this means that the background is changing while the tensor wave amplitude not. 

In order to study the cosmological stochastic background, the operator (l42t can be specified for a RW fiat model, 
i.e., K = in Eg. (ITT]) , and then the propagation equation dTH)) becomes 

h+[3H + 20] h + k^a-^ h = , (43) 

being □ = df + 3Hdt and k — kik^ the comoving wave number. It is worth stressing that Ea. (IT51) applies to any J{R) 
theory whose conformal transformation can be defined as e^"^ = f'{R) and, thus, the solution, i.e., the GW amplitude, 
depends on the specific cosmological background (i.e., a{t)) and on the specific theory of gravity {i.e., 0(i)). 

We start by studying a non-analytical /(i?) model with a power-law scale factor. We assume a power-law behavior 
for the scale factor, i.e., a = ag [t/to]^ (we restrict ourselves to the case x < 1), and the non-analytical form of f{R) 
given in Eq.dl]). Using the relation (j)(t) = (1/2) In/' and Rrw = 6a;(2x — l)/t^ (for K = 0), we get 

0(i) = 1 In [l + = i In 1 + nq{6x{2x - 1))""^ t-2("-i) . (44) 



It is straightforward to see that during the radiation-dominated era {x ~ 1/2), 0(i) vanishes identically. Thus Eq. p3)) 
reduces to its GR form and the standard result for the propagation of GWs follows, i.e., h on sm(kri) / krj , where rj is 
the conformal time, defined by the relation dt = adrj. 

For X ^ 1/2, the time derivative of the field (j) is given by: 



n(l-n)g[6a;(2x- 1)]" t^-^" _ (n ^ 1) 1 

1 + nq[6x{2x - 1)]"-' i-2(„-i) " t (t/t,)2("-i) - 1 ' 
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where we defined = {n\q\)^^'^^"'~^'' |6a;(2x — l)p. We remark that the second equality is vaUd for n in the form 
{2i + 1)/(2to + 1), as required by the model. Thus, it is easy to see that for t <C t*, — (1 — n)/t. Since iJ ~ 1/t, in 
this regime (p and H are comparable and we expect the f{R) action to affect the propagation of GWs. When t ^ t^., (f) 
becomes very large (and eventually diverges for t — t^) and the corrections to the EH action dominate the dynamics. 
Finally, when t ^ t^,, we have that |(/)| ^ and the standard behavior is recovered. 

Let us now briefly discuss the GW propagation in the regime t <C <*, when = (1 — n)/t. In this case, Eg. pSl) can 
be recast in the form 

h P 

h+[3x + 2{l-n)]- + —h^O . (46) 
t 

It is convenient to switch to the time variable u = krj. Eq. (1461) thus rewrites as 

(Ph 2(x-n + l)dh ^^^^ 



du'^ (1 — x)u du 

The main advantage of this variable choice is that the explicit dependence on the wave number k disappears from 
the dynamical equation. Furthermore, since rj represents the cosmological horizon, it is easily understood that u = 1 
roughly corresponds to the time of horizon crossing, and that the regime u <^ 1 (m 3> 1) represents a wave that is far 
outside (inside) the horizon. The general solution to this equation is 

Hu) ^ [A,Uu) + A.Y^iu)] , « = ^77^ ^. (48) 

2(1 — x) 

where and Ya are the Bessel functions of the first and second kind, respectively, and the integration constants Ai , 
A2 are fixed by the initial conditions. A useful insight can be gained analyzing the asymptotic behavior of the above 
solution for it ^ 1 and u ^ 1. In the range of interest (2 < n < 3 and < a; < 1), we have that a > 3/2 and we can 
safely use the following expansions for the Bessel functions: 

/ 2 / air 7r\ / 2 / air 7t\ 

Ja[u)~\ — cos u ;r " T ' Ya{u) \ — sm u tt ^ T I ' (49b) 

VTTuV 2 4/ \ nu \ 2 4/ 

Thus the two independent solutions in Eq.(j3S]) are both well-behaved for u — 0; in particular, the solution Jq,(m) — 
0, while u^Yalu) tends to a constant value. The first derivatives of both solutions are also vanishing in u = 0, as 
it can be shown by using the relations between the derivative of Jq in terms of Ja-i and Ja+i (and similarly for 
Ya). On the other hand, the behavior of the solution when the wave is well inside the horizon (u ^ 1) shows some 
interesting features. In fact, in this regime both solutions oscillate with an amplitude that grows with time like 
(we recall that a > 3/2); this is in striking contrast with the usual behavior that is found in standard GR, where the 
amplitude decays with time {e.g., h oc siii{u)/u during the radiation-dominated era). This is related to the fact that 
the coefficient of the first-derivative term in Ea. (l47|) . that is positive in standard GR, becomes negative for 71 > a; + 1, 
or, equivalently, for a > 1/2. 

As already said, during the radiation-dominated era, vanishes identically and the analysis above does not apply. 
In principle, however, the propagation of GWs during the matter-dominated era (a; — 2/3) could be modified with 
respect to the standard picture, depending on the value of i*, that we are now going to estimate. As discussed in the 
previous section, the value of \q\, or equivalently of Lq, is bounded from above by the requirement of the existence 
of the matter-dominated era itself. The exact value of L^^"^ depends on n, but only weakly, so that we can safely 
consider that Lg < 70 pc. Then we find that < 400 yr; since the matter-dominated era began at t ~ 10^ yr, we can 
conclude that during this epoch the contribution of the <p term is negligible and thus the propagation of GWs is 
standard. 

To conclude, we briefly turn our attention to the GW propagation during a de Sitter phase. We then hypothesize 
an exponential behavior for the scale factor of the Universe a = d e^^*, with Hj > 0, as done in Section ITVl In this 
case, using the relation (f){t) = (1/2) In/', one can argue that 

= i In(ng) + 1 (n - 1) In [l2H] + GKe'^"'* /d"^] , (50) 

and, of course, if we focus on the case of vanishing spatial curvature K — 0, one gets (j) — 0. Since, as it can be seen 
from Eq. p5|) . the modifications to the GR action only enter the GW dynamics through a term proportional to (p — 0, 
this implies that their propagation during a de Sitter phase is the same as in standard GR. 
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VI. CONCLUDING REMARKS 

In this work, we have considered a modified theory of gravity of the form f{R) — R + qK"', where both q < and 
2 < ri < 3 are free parameters. 

The present analysis has demonstrated how, in contrast with analogous treatments [sj, the non-analytic case in 
exam is not excluded by Solar-System tests. Instead, using the measurements of the Earth period, a lower bound for 
the characteristic length scale Lq = |g|i/(2n-2) ^j^g proposed model can be found, shedding light on the viability 
of the theory in the weak- field limit. Such constraint confines non-Newtonian effects to be relevant outside the Solar 
System only, unless n is fine-tuned to be very close to 2. We underline that our proposal for a suitable value of 
the characteristic length could induce, in principle, relevant modifications in larger-scale astrophysical systems, e.g., 
stellar clusters. 

The model has also been implemented in a cosmological framework, showing the possibility to recover the main 
phases of the Universe evolution, i.e., the radiation- and matter-dominated era. We have assumed a power-law 
behavior of the cosmological scale factor and have shown how these two phases are reproduced in the relevant limits. 
In fact, the radiation-dominated evolution is recovered (irregardless of the spatial curvature) in the asymptotic limit 
toward the singularity. In the same limit, a phase of power-law inflation was also found for vanishing curvature. On 
the other hand, in the limit of large times, the matter-dominated phase is reproduced in the case of zero curvature. 
As a main result, an upper bound of ^ 70 pc for the characteristic length was found by guaranteeing the existence 
of such an era of the evolution of the Universe. The presence of this cosmological bound and of the Solar-System 
constraints, has allowed us to define a restricted region of the parameter space in which the model can be implemented 
both in the weak-field limit and in the cosmological framework. 

Moreover, we have shown how an exponential early evolution of the Universe is still present in the extended f{R) 
scheme. The expansion rate was found to be much smaller than the standard value by several orders of magnitude due 
to the presence of a minimum length scale for the model. This allows to have the inflationary scale arbitrarily close 
to the Planck one without spoiling the current limits on the tensor-to-scalar ratio r, since the latter is proportional 
to the expansion rate and thus gets correspondingly suppressed. More generally, we have found that a tiny value of 
the tensor-to-scalar ratio is expected in the extended framework unless the exponent n is fine tuned to be unnaturally 
close to 2 (but still different from 2). The observational implication is that, if the gravitational action is of the form 
considered here, tensor modes will not be detected even by next-generation CMB-polarization experiments. On the 
other hand, a detection of tensor modes in the near future would falsify the present model. 

The issue of the GW propagation on a RW background has also been addressed, assuming a power-law as well as 
an exponential behavior of the cosmological scale factor. We have shown the presence of a typical time scale over 
which /(i?) corrections can be relevant, which implies that the term due to the generalized framework can always 
be neglected during the matter-dominated era. On the other hand, the corrective term vanishes identically both in 
the radiation-dominated era and in the de Sitter phase. This way, the GW evolution results to be, for all practical 
purposes, the same as in standard GR. Combined with the suppression of primordial tensor modes, this implies that 
the inflationary GW background will be below the detection threshold of present and future interferometers, and thus 
represents another observational test for the model. On the contrary, different production mechanisms {e.g., cosmic 
strings) could still, in principle, yield a cosmological GW signal of sizable intensity. 
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